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In contrast to the nuclear structure along the beta-stability line which has been well 
studied both experimentally and theoretically, the yet unknown structure of drip-line nuclei 
is currently of great interest From the theoretical point of view, spectroscopy of 

exotic nuclei offers a unique test of those components of effective interactions that depend 
on isospin degrees of freedom. Since the parameters of interactions used in the usual mean- 
field calculations are determined so as to reproduce the properties of beta-stable nuclei, these 
parameters may not always be optimal around particle drip lines due to (often dramatic) 
extrapolations involved. 

The nuclei discussed in this study are the Sulfur isotopes, especially the neutron-rich ones 
with N~28. This particular choice was motivated by recent experimental and astrophysical 
interest in this mass region From the perspective of the spherical shell model, the 

underlying proton wave functions are described in terms of the sd shell-model space, while 
the main neutron components originate from the fp shell-model space. Such a schematic 
classification, however, easily breaks down due to the strong core polarization effect (i.e., 
the appearance of static shape deformations associated with the core-breaking excitations). 
Experimentally, deformed states in magic nuclei are known in many cases (see the recent 
review [|^). Those intruders sometimes appear very low in energy and, in a few cases, they 
become ground states. 

In the sd region, the neutron-rich nuclei with N~20 are spectacular examples of coexis- 
tence between spherical and deformed configurations. A classic example of a magic nucleus 
with a deformed ground state is 20^85 which has a very low-lying 2+ state at 886 keV 
PI and an anomalously high value of S2n- Calculations based on the deformed mean- field 
theory predict deformed ground states around '^^Mg and explain them in terms of neutron 
excitations to the I/7/2 shell A similar conclusion has been drawn in the shell model 

calculations in the {sd){fp) model space [0,0 and in the schematic analysis of ref. 

Another region of unexpected collectivity are the I/7/2 systems around 24Cr24. Naively, 
the main features of these nuclei should be well reproduced in terms of a single-j shell-model 
picture governed by pairing interaction. However, many properties of the observed states 
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cannot be accounted for by the results of the empirical (I/7/2)" shell model calculations |T4 



and the extension to the full (/, p) configuration space is necessary [|T^,|I^. Several nuclei 
around ^'^Cr exhibit rotation-like level spacings |T7| , p!8[ and the self-consistent calculations 
|T^-^ yield deformed ground states. In particular, ^'^Cr is calculated to be prolate-deformed 
with quadrupole deformation /?2=0.28 

Experimentally, little is known about the neutron-rich nuclei around '^^S. Recently, (3- 
decay properties of ^^S and 45-47qj have been studied in GANIL l2^,§]. The half-life of ^^S was 



found to be Ti/2=123±10 ms. (Observation of ^^Si, ^^'^^p^ 48g^ ^^^^ 51qj j^^^g j^gg^ reported 



in ref. p3[.) The structure of exotic neutron-rich nuclei with 10<Z<20 will soon be studied 
at GSI using the fragmentation reaction ^^Ca on ^Be at relativistic energies (f/c~0.5) 

Theoretical information on the light N^28 nuclei is also very scarce. The stability of 



highly neutron-rich Si isotopes was investigated in ref. in the spherical HF framework 
with various Skyrme interactions. The results of recent large-scale mass calculations using 
the finite-range droplet model (FRDM) |^ or the extended Thomas-Fermi with Strutinski- 
integral model (ETFSI) suggest the presence of deformation in some N=28 isotones 
(see, e.g.. Table IV of ref. ^). Based on these calculations and on the theoretical analysis of 
/3-decay properties of '^''S, it has been suggested in ref. |^ that the influence of the spherical 
shell N=28 is weakened below '^^Ca. The microscopic structure of neutron-rich N~28 nuclei 
is important in the astrophysical context. Indeed, the neutron-rich N^28 nuclei are crucial 
for the nucleosynthesis of the heavy Ca-Ti-Cr isotopes 

To shed some light on the physics of exotic neutron-rich nuclei with N:^28, we performed 
calculations based on the self-consistent mean-field theory, namely the Skyrme-HF model 
and the relativistic mean-field (RMF) model. In this Letter we report results for the Sulfur 
isotopes. The results for the Silicon and Argon isotopes will be published in a forthcoming, 
more detailed paper. 

We perform the Skyrme-HF calculations by discretizing the energy functional on a three- 
dimensional Cartesian spline collocation lattice, which provides a highly accurate alternative 
to the finite-difference method The structure of the resulting lattice representation is 



suited for vector and parallel supercomputers, and the method allows for highly modular 
programming where the order of the splines can be defined as an input parameter. Equations 
of motion are obtained via the variation of the lattice representations of the constants of 
motion, such as the total energy. It is worth noting that no self-consistent symmetry has 
been imposed in the calculations. The details of our method have been published in Ref. pS 



In this work we use the Skyrme interaction SIII The exchange part of the Coulomb 
interaction is taken to be in the Slater form. In addition, we use the simple scaling of 
the nuclear mass to approximately correct for the center-of-mass motion. The calculations 
were performed in the cube of the size (20 fm)^. In all the cases considered, the resulting 
HF minima turned out to correspond to reflection-symmetric shapes with three symmetry 
planes, i.e., (x) = {y) = {z) = Q and {xy) = {xz) = {yz) = 0. 

The basic building blocks in the relativistic mean-field approach ||3^ are the baryons 
(protons and neutrons) and the a—, u—, and p— mesons. The cr-meson is assumed to move 
in a non-linear potential 

U (a) = ^m^a^ + ^g^a^ + ^g^a^ (1) 

In the present work we have employed the recent set of Lagrangian density parameters NL- 
SH ref. This set has been claimed to be particularly successful in describing properties 
of very neutron-rich systems. The relativistic equations of motion are derived by means of 
the variational principle. The resulting Dirac equation for the baryons and the Klein-Gordon 
equations for the mesons are solved using the basis expansion method. In this method the 
small and the large components of the Dirac spinor and the meson fields are expanded in 
terms of the axially symmetric-stretched harmonic oscillator basis with oscillator frequency 
^u;=41A~^/'^MeV. (twelve deformed oscillator shells for neutrons and protons were used). 
In the RMF calculations, only reflection-symmetric axial shapes were considered. 

It is known |^ that in drip-line nuclei, pairing interaction leads to the scattering of 
nucleonic pairs from bound states to continuum. In the mean-field-|-BCS model this leads 
to the presence of unphysical "particle gas" surrounding the nucleus. The deficiencies of the 



standard treatment of pairing around drip lines can be cured by means of the Hartree-Fock- 
Bogolyubov (HFB) method with a reahstic pairing interaction in which the wave functions 
of occupied quasiparticle states have correct asymptotic behavior Unfortunately, the 

three-dimensional HFB code is not available at present. Consequently, in our study we used 
the "constant-gap" approximation with a strongly reduced pairing, A„=Ap=75 keV and 
200 keV in the RMF and HF models, respectively. We also performed the RMF calculations 
with A„=Ap=500 keV. The results obtained were very similar to those with A„=Ap=75 
keV. 

The ground-state minima of the Sulfur isotopes, in most cases, can be associated 
with deformed intrinsic states. In order to compare various variants of calculations and 
to relate them to previous work, the standard quadrupole deformation parameters (32 



and 7 were extracted. Firstly, the two quadrupole moments, g2o = y 167r/5(r^F2o) and 
q22 = y'87r/15(r^ {Y22 + Y2-2)) were expressed in terms of the polar coordinates Qo and 7 
3|,|35 

q2Q = Q0COS7, q22 = --^QoSin7. (2) 

Using Eq. (g), the 7-value can be determined. The quadrupole moment of proton distribu- 
tion, Qo, can be written in terms of /^f by means of relation 

Q'o = Z(rJ)/3f. (3) 

Similarly, one can extract quadrupole deformations of neutron {/32) and mass (/?^) distri- 
butions. Since, in most cases, the equilibrium shapes predicted by the HF model are axial, 
we adopted the standard convention (i.e., /32>0 (prolate) for 7~0° and /?2<0 (oblate) for 
7~60°). 

The two-neutron separation energies for the Sulfur isotopes, S2n, are displayed in Fig. |1| 
as a function of neutron number. Our HF and RMF results are compared with the predic- 
tions of the FRDM and ETFSI models, and experimental data. In general, the agreement 
between the models themselves and between theory and experiment is good. In the HF 



model the nucleus ^^S (N=36) is two-neutron-unstable. The RMF calculations yield more 
neutron binding; the isotopes ^^'^^S are predicted to lie inside the two-neutron drip-line. 
The systematic model agreement for S2n should not imply that the intrinsic structures of 
neutron-rich Sulfur isotopes are also similar in all models presented in Fig. |^. As discussed 
below, this is not the case. 

The calculated mass quadrupole deformations of even-even Sulfur isotopes are shown in 
Fig. 1^ as a function of neutron number. (Usually, calculations yield more than one energy 
minimum. In such situations, deformations and excitation energies of excited states are also 
displayed.) The deformation pattern for even-even isotopes ^s-ssg jg fairly similar in the HF 
and RMF models: prolate ground states in 28-32,38g^ oblate minimum in ^"^S, and spherical 
shape in the magic N=20 nucleus '^^S|. However, the differences show up for the heavier 
isotopes. 

The nuclei ^^'^^g ^.^^ predicted by the RMF model to have prolate ground states with 
/92'^0.25; the oblate minima with /32~-0.16 lie about 4 MeV higher in energy. In the HF 
calculations, prolate (/52~0.25) and oblate (/52~-0.24) configurations are practically degen- 
erate. The FRDM and ETFSI models give prolate-deformed ground states with (/52~0.25), 
in agreement with our findings. 

There is no consensus regarding the equilibrium shape of the N=28 nucleus ^^S. According 
to RMF calculations, ^^S has a well-deformed prolate ground state with /32=0.31. The oblate 
minimum with /?2=0.16 appears to lie 0.8 MeV higher. The spherical saddle point lies even 
higher at £'*=2.8 MeV. On the other hand, according to the HF model, ^^S is a gamma- 
soft system with a small quadrupole deformation /52~0.13. The very shallow, well-deformed 
prolate minimum (/?2=0.28) analogous to the RMF ground state lies at £'*~1.5 MeV. The 

•^It is to be noted here that RMF predicts a slightly deformed shape for ^^S. This can be attributed 
to very weak pairing in the calculations. If a stronger proton pairing gap is assumed, Ap=l MeV, 
the equilibrium shape becomes spherical. 
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FRDM and ETFSI models give spherical and oblate {^2=^0.26) ground states, respectively. 

The heavier Sulfur isotopes with N>28 are consistently calculated to be prolate-deformed 
in the HF, RMF, and FRDM models. The oblate minima are higher in energy, but with 
increasing neutron number, the prolate-oblate energy difference is reduced. The RMF ap- 
proach yields systematically larger quadrupole deformations compared to HF. For instance, 
according to the RMF predictions, ^^S is a well-deformed system with /?2=0.22. According 
to the HF model, it is a transitional nucleus with /?2=0.13, 7=10°. Interestingly, the ETFSI 
approach favors a strongly deformed oblate configuration with /?2=^0.25. 

For the N=28 isotone, ^^Si calculations yield well-deformed oblate ground states with 
p2 ranging from -0.32 (FRDM) to -0.22 (HF). Also, for the Sihcon isotopes, the RMF 
model gives slightly more neutron binding than HF. According to RMF, the nucleus ^'^Si is 
still inside the two-neutron drip-line, while it is unstable in the HF model. Another N=28 
isotone, "^"^Ar, is predicted to be spherical (RMF, FRDM) or oblate (HF) with /?2=-0.14. 
The RMF model gives an oblate minimum with /52=-0.10 at very low excitation energy. 
This indicates the transitional character of '^^Ar. 

Figure |^ displays the representative single-particle neutron levels as functions of 
quadrupole deformation ^2- The Nilsson diagram was obtained using a deformed Woods- 
Saxon model with a Chepurnov set of parameters the resulting single-particle energies 



are close to those from the HF+SIII model. (In the RMF-NL-SH calculations the order of 
Si/2 and (i3/2 shells is reversed.) The neutron and proton single-particle levels in the ground- 
state of ^^S obtained in HF and RMF models are also displayed in Fig. ^. The deformed 
shape in ^^S results from a subtle interplay between the deformed gaps at Z=16 and N=28, 
and the spherical N=28 gap. In the RMF model, the spherical N=28 gap is completely bro- 
ken; the neutron intruder orbital [321]l/2, originating from the 2^3/2,1/2; I/5/2 shells above 
the N=28 gap is occupied. In the HF calculations the energy distance between the deformed 
[321]l/2 and [303]7/2 levels is only 1.6 MeV. 

The predicted root-mean-square (rms) neutron and charge radii of even-even Sulfur iso- 
topes are illustrated in Fig. ^. The results of HF and RMF models are fairly similar. There 



appears a systematic shift of ~0.08 fm between the HF and RMF results for the rms charge 
radii; the experimental data for 32,34,36g jjg between. As expected, due to skin effect, 
rms radii increase when approaching particle drip-lines. (Since, in our study, the pairing 
correlations are practically neglected, the undesired "particle gas" effect mentioned above is 
not present.) 

The difference between neutron and proton deformations, Aj32=l32 — ^2, is illustrated 
in Fig. |[ It is seen that when approaching the neutron drip-line, the values of 
systematically smaller than those of the proton distribution. An opposite effect is seen 
around the proton drip-line. The largest difference, |A/52|~0.10, is obtained in the RMF 
model for ^^S. (As discussed above, the deviation for ^^S disappears if the stronger proton 
pairing gap is assumed.) The behavior of A/32 can be partly attributed to the isotonic 
behavior of rms radii in Fig. H (see the definition (^ of ^2)- Indeed, the value of Qo depends 
both on the angular anisotropy and the radial dependence of the nucleonic density. In the 
drip line nuclei, due to spatially extended wave functions, the "radial" contribution to Qo 
might be as important as the "angular" part. Moreover, it can strongly depend on particle 
number. The insert in Fig. ^ illustrates the dependence of the Q^/Qo ratio on the N/Z ratio. 
Up to neutron number N=26 the HF and RMF results are very similar. The dips in Q^/Qo 
in the HF calculations appear at shell and subshell closures (i.e., at N=14, 20, 28, and 32). 
At these particle numbers the neutron distribution has a tendency to be more spherical. 
The lack of shell fluctuations in Q^/Qo in the RMF model at N>26 is consistent with 
the calculated large prolate deformations. Interestingly, above N=28, A/?2(RMF) decreases 
steadily with neutron number, while the Q"/Q^(RMF) ratio intersects the N/Z line (rigid 
geometric limit) only at N>34. 

In summary, our calculations suggest strong deformation effects in the region around 
"^^S due to the l/7/2~^/p core breaking. Large differences between neutron and proton 
quadrupole moments are obtained in the RMF approach in deformed Sulfur nuclei far from 
stability. Such differences might have interesting consequences for the quadrupole isovector 
modes in drip line systems. According to the RMF+NL-SH model, the N=28 nuclei ^^S 
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and Si are well deformed in their ground states. On the other hand, the HF+SIII model 
predicts smaller deformations, and ^^S is calculated to be a soft transitional system. The 
new experimental data around "^^S will certainly be very helpful in pinning down the question 
of quadrupole collectivity in this mass region. 
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FIGURES 

FIG. 1. Two-neutron separation energies of the even-even Sulfur isotopes calculated with the 



HF and RMF models. They are compared with the results of the FRDM |2|] and ETFSI ^ 
models, and experimental data. 

FIG. 2. Quadrupole mass deformations of the even-even Sulfur isotopes calculated with the 
HF (top) and RMF (bottom) models. The dots connected by a solid line correspond to ground-state 
deformations. The empty circles indicate excited configurations (with excitation energies given in 
MeV). 

FIG. 3. Left: Woods-Saxon single-particle neutron levels as functions of 02- The orbitals are 
labeled by means of 0, and vr {7r=+, solid line; 7r=-, dashed line) quantum numbers. Right: The 
shell structure of ^^S as predicted by the HF and RMF models. 

FIG. 4. Root-mean-square radii of the neutron distribution (top) and charge distribution (bot- 
tom) of the even-even Sulfur isotopes calculated with the HF and RMF models. The radii are 
shown relative to the average (liquid-drop) value of y/3/5Ro, Ro = 1.2A^/^ fm. The experimental 
data for 32,34,36 g ^^^^ taken from ref. [pGl 



FIG. 5. Difference A/?2=/32 ~ even-even Sulfur isotopes calculated with the HF and 

RMF models as a function of neutron number. The insert shows the ratio of neutron and proton 
quadrupole moments versus N/Z. The large-deformation excited state in ^^S calculated in the HF 
model is indicated by means of an "(•)" symbol. 
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